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COMPUTING  EQUILIBRIA  VIA  NONCONVEX  PROGRAMMING 

by 

Jonathan  F.  Bard 
James  E.  Falk 


1.0  Introduction 

In  this  paper  we  investigate  a procedure  for  computing  equilibria 
from  the  vantage  point  of  mathematical  programming.  A competitive  model  of 
an  economy  will  serve  as  the  basis  for  the  discussion  although  a variety  of 
contexts  would  have  been  equally  suitable.  Other  types  of  equilibrium  prob- 
lems, such  as  those  arising  in  traffic  network  analysis,  have  direct  con- 
ceptual and  analytic  counterparts  to  those  found  in  economics,  and  are 
hence  amenable  to  the  same  solution  techniques. 

A state  of  equilibrium  exists  when  competing  or  opposing  forces  are 
brought  into  balance.  One  of  the  major  themes  of  economic  theory  is  that 
the  behavior  of  a complex  economic  system  can  be  viewed  as  an  equilibrium 
arising  from  the  interaction  of  a number  of  economic  units,  each  motivated 
by  their  own  special  interests.  General  equilibrium  theory  ([2],  [22], 

[24])  seeks  to  determine  the  point  at  which  this  balance  can  be  struck,  and 
in  so  doing  focuses  on  the  interrelationships  that  exist  among  the  markets 
for  goods  and  services  in  the  economy.  The  analysis,  however,  is  carried 
out  in  terms  of  individual  decision  makers  and  commodities  rather  than  in 
terms  of  aggregates.  The  fundamental  questions  that  general  equilibrium 
theory  attempts  to  answer  are  the  same  as  those  posed  in  macroeconomic 
theory:  given  different  economic  environments,  what  goods  will  the  economy 
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produce,  how  will  these  be  produced,  and  who  will  obtain  them?  But,  where 
macroeconomics  provides  answers  in  terms  of  aggregates,  general  equilibrium 
theory  provides  answers  in  terms  of  the  individual  consumers,  producers, 
and  commodities  making  up  these  aggregates. 

Consider,  for  the  moment,  a model  in  which  m consumers  are  engaged 
. in  the  exchange  of  commodities  which  they  initially  own  and  in  which  production 

or  supply  is  ignored.  Suppose  there  are  n goods  in  the  economy  and  that 
» each  of  the  consumer's  preferences  is  represented  by  a utility  function.  A 

bundle  of  goods  x is  preferred  to  a bundle  x'  by  consumer  i (i  = l,..,,m) 
if  and  only  if  u^(x)  > u^Cx')  where  the  utility  function  u^  : Rn  •>  R is 

generally  assumed  to  be  strictly  concave  and  continuous.  Let  p e R°  be  the 
vector  of  prices  for  the  n goods.  The  demands  of  the  ith  consumer  are 
determined  by  the  solution  to  the  following  problem: 

maximize  u.(x) 

1 

subject  to  <p,x>  < <p,wX> 
x > 0 

where  w1  is  the  initial  wealth  or  resource  endowment  of  the  ith  consumer, 
i = l,...,m  . We  shall  assume  that  the  solution  vector  for  this  problem, 
d^-(p)  , can  be  written  as  a continuous  function  of  the  prices  p . The 
individual  trader's  excess  demand  function  is  d^p)  - w*  (i  = l,...,m) 

• and  will  be  denoted  by  g^(p)  . The  excess  demand  will  be  positive  for 

those  commodities  whose  stock  he  wishes  to  increase  by  exchange  and  negative 
for  the  remaining  items.  If  it  is  assumed  that  all  purchases  are  to  be 
financed  solely  by  the  sale  of  assets,  then  individual  budgetary  constraints 
lead  to  the  following  identity: 

I 

P-id^(p)  + •••  + p d^(p)  = p w^  + •••  + p w1  . (1) 

\ j-j.  rinii  nn 

The  market  excess  demand  function  g : Rn  ->  R°  is  simply  the  sum  of  the 
individual  excess  demand  functions 


m i i 

g(p)  = E (d  (p)  - w1)  . 

i=l 


An  equilibrium  price  vector  p*  is  one  for  which  all  of  the  market 
excess  demands  are  less  than  or  equal  to  zero  with  a zero  price  for  any 
commodity  whose  excess  demand  is  strictly  less  than  zero.  This  leads  to 
the  formulation  of  the  complementarity  problem  ([8],  [17]): 

g(p)  < 0 , (2) 

P > 0 , (3) 

<P»g(p)>  = 0 , (4) 

whose  solution  p*  will  be  the  focal  point  of  this  paper.  Condition  (4), 
known  as  the  Walras  Law  [27],  is  obtained  by  aggregating  individual  excess 
demands  to  create  the  set  of  market  excess  demand  functions,  and  then  put 

in  the  form  of  (1).  Walras'  law  holds  for  all  price  vectors  p whether 

they  be  in  equilibrium  or  not. 

We  note  here  that  production  may  easily  be  incorporated  in  this 
model  by  either  replacing  or  augmenting  the  ith  consumer's  initial  wealth 
w1  by  a supply  function.  For  individual  i , this  function  relates  the 
prevailing  market  prices  p to  the  quantity  of  goods  produced. 

A number  of  persons,  including  Nash  [20],  Arrow  and  Debreu  [1],  and 
Kuhn  [16]  have  studied  the  existence  problem  of  the  competitive  model  from 
the  standpoint  of  combinatorial  topology.  The  first  algorithms,  however, 
actually  designed  for  computing  economic  equilibria  were  developed  by  Scarf 
[23],  and  were  based  on  a procedure  for  approximating  a fixed  point  of  a 
continuous  mapping.  More  recently,  Wilson  [28]  and  Elken  [10]  have  exploited 
path  methods  in  the  pursuit  of  greater  computational  efficiency.  In  a 
slightly  different  vein,  Lemke  [17]  offered  some  constructive  proofs  relating 
to  the  existence  of  equilibrium  points  for  bimatrix  games.  His  work  strongly 
suggested  a computational  scheme  for  models  with  linear  excess  demand 
functions. 

This  paper  presents  an  alternative  procedure  for  computing  equilibria 
for  a class  of  problems  where  the  excess  demand  function  or  its  logical 
equivalent  has  an  explicit  representation  that  can  be  converted  to  a separable 
form.  Solutions  are  obtained  by  first  recasting  the  complementarity  problem 
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into  a nonconvex  minimization  problem  whose  optimal  value  or  best  upper 
bound  is  known  at  the  outset,  and  then  using  Falk's  [12]  algorithm  to  locate 
a global  solution.  This  allows  us  to  go  beyond  the  common  linear  formula- 
tions of  an  economy  or  network  (e.g.,  see  Eaves  [8],  Negishi  [21],  or  Asmuth, 
Eaves,  and  Peterson  [4])  which,  in  spite  of  their  outward  simplicity,  must 
appeal  to  rather  complicated  algorithms  if  solutions  are  to  be  obtained. 


The  algorithm  which  we  subsequently  describe  and  use  as  an  alternative 
for  solving  (2)  - (4)  is  based  on  a branch  and  bound  philosophy,  and  as  such, 
computes  a convergent  sequence  of  upper  and  lower  bounds  on  the  optimal  value 
of  the  problem.  In  our  case,  however,  because  the  best  upper  bound  on  the 
objective  function  is  known  to  be  zero,  the  amount  of  work  necessary  to 
achieve  convergence  is  significantly  reduced.  The  usual  requirement  of  find- 
ing a point  that  yields  equality  between  the  best  upper  and  best  lower  bounds 
is  replaced  by  the  simpler  requirement  of  finding  any  point  that  yields  an 
objective  value  of  zero. 

In  the  next  section,  the  complementarity  problem  is  reformulated  as 
a nonconvex  minimization  problem  whose  solution  yields  the  desired  equilibrium 
vector.  Next,  the  method  is  applied  to  a number  of  sample  problems  and  our 
computational  experience  is  detailed.  Here  we  see  that  the  results  are 
obtained  in  a surprisingly  small  number  of  iterations  of  the  algorithm. 

2 Reformulation  of  the  Complementarity  Problem 

In  the  complementarity  problem  derived  above,  there  is  no  objective 
function  to  be  optimized.  Indeed,  in  many  complex  economic  equilibrium 
problems  there  does  not  appear  to  be  a "natural"  objective  function  whose 
optimization  yields  prices  and  quantities  in  equilibrium  (see,  e.g..  Scarf 
[24]). 

In  spite  of  this,  consider  the  following  "artificial"  minimization 
problem  P: 

v*  = min{-<p,g(p)>  : g(p)  < 0 , p > 0}  (P) 

Now  let  p*  be  a solution  of  the  complementarity  problem  (i.e.,  p*  is  a 
vector  of  equilibrium  prices).  Then  p*  is  feasible  to  problem  P,  and 


r- 
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yields  a value  of  0 to  the  objective  function.  Since  this  objective 
function  is  greater  than  or  equal  to  zero  at  all  feasible  points,  v*  = 0 . 
Conversely,  it  is  clear  that  any  solution  of  (P)  for  which  v*  = 0 must 
be  a vector  of  equilibrium  prices. 


Problem  P is  of  a nonconvex  nature,  and  in  general,  no  suitable 
technique  exists  for  the  determination  of  a global,  rather  than  a local 
solution;  however,  if  each  excess  demand  function  g , i = l,2,...n  , is 
separable,  i.e.. 


81(P) 


I g,,(p.)  1 = 1>2 


j=l 


ij 


,n 


and  each  g_  is  continuous,  then  (P)  can  be  written  as  a separable  pro- 
gramming problem  whose  approximate  global  solution  can  be  obtained  with 
arbitrary  precision. 


We  now  formulate  an  equivalent  problem  with  a different  objective 
function  but  the  same  constraint  region  whose  optimal  value  is  equal  to 
that  of  (P).  The  equivalent  problem  P'  is 

min{£  min(pi,-gi(p))  : p > 0 , g(p)  < 0}  (Pf) 

i 

Rewriting  the  objective  function  in  (P1),  we  get  the  desired  result: 

min  {£(min(0,-g  (p)-p  ) +p.)} 
g(p)^0  i 
p>0 


= min{£(min(0,w.)  + p.)} 
i 

g(p)^o 

w+p+g(p)  = 0 

p>0 


= min{ ^ (min(0 jW^)  + p^)} 


(S) 


J8ij(pj>  i 0 

w.+p .+Fg. . (p . ) = o 
i 

p>0 


i = 1,2, . . . ,n 


where  the  w ' s will  be  referred  to  as  auxiliary  variables. 
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Problem  S is  still  a nonconvex  programming  problem,  but  its  separable 
structure,  created  at  the  expense  of  a twofold  increase  in  dimensionality, 
makes  its  mathematics  much  more  tractable.  The  traditional  method  for  treat- 
ing separable  problems  involves  calculating  piecewise  linear  approximations 
of  the  associated  functions  and  applying  a modification  of  the  simplex  method 
to  the  resulting  problem  (see,  e.g.,  Miller  [19]).  The  modification  amounts 
to  a restriction  on  the  usual  manner  of  selecting  variables  to  exchange  roles 
(basic  to  nonbasic  and  vice  versa)  and  will  yield  a local  but  not  necessarily 
a global  solution  of  the  approximating  problem. 

An  algorithm  for  finding  global  solutions  of  nonconvex  separable 
problems  was  developed  by  Falk  and  Soland  [13]  and  Soland  [25].  The  method 
is  based  on  the  branch  and  bound  philosophy  and  yields  a (generally  infinite) 
sequence  of  points  whose  cluster  points  are  global  solutions  of  the  problem. 
The  implementation  of  the  method  is  limited  by  the  necessity  of  computing 
convex  envelopes  [11]  of  the  functions  involved  although  a number  of  applica- 
tions have  been  shown  possible  when  these  functions  exhibit  special  structures 
(e.g.,  concave  or  piecewise  linear). 

The  inherent  limitations  that  special  problem  structures  impose  have 
been  overcome  by  the  introduction  of  two  algorithms  independently  developed 
by  Beale  and  Tomlin  [5]  and  Falk  [12].  For  this  paper , we  have  used  the 
programming  code  MOGG  based  on  the  algorithm  proposed  by  Falk  and  written  by 
Grotte  [15],  to  solve  a number  of  equilibrium  problems.  The  results  are 
presented  in  the  next  section. 

3.0  Computational  Experience 

A variety  of  equilibrium  problems  have  been  studied  to  test  the 
approach  outlined  above.  The  first  is  a multicommodity,  transshipment  prob- 
lem defined  on  an  affine  network,  taken  from  Asmuth,  Eaves,  and  Peterson 
[4]  who  used  Lemke's  algorithm  [17]  to  obtain  a solution.  The  second  involves 
a simple  competitive  market  comprising  three  producers,  three  consumers,  and 
three  commodities.  The  supply  and  demand  functions  in  this  economy  are  given 
a piecewise  linear  formulation,  and  three  equilibrium  points  are  known  to 
exist.  The  third  problem  is  identical  to  the  second  except  a majority  of 
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the  piecewise  linear  functions  have  been  recast  as  continuous,  smooth 
functions.  The  fourth  problem  provides  an  example  outside  the  context  of 
economics,  and  is  derived  from  a 3 node,  4 arc  traffic  network  whose  equiva- 
lent excess  demand  function  is  both  nonlinear  and  nonseparable . 

The  algorithm  itself  is  based  on  branch  and  bound  techniques  and 
considers  subsets  of  a linear  polyhedron  containing  the  feasible  region  of 
(S).  A lower  bound  on  the  optimal  value  of  the  problem  is  found  by  minimizing 
the  objective  function  over  each  of  these  subsets  and  selecting  the  smallest 
value  obtained.  A check  for  the  solution  is  made  which,  if  successful, 
yields  a global  solution  of  the  piecewise  linear  approximation  to  (S).  If 
the  check  fails,  the  subset  corresponding  to  the  smallest  lower  bound  is 
further  subdivided  into  either  two  or  three  new  linear  polyhedra  and  the  pro- 
cess continues  as  before  with  new  and  sharper  bounds  being  determined.  The 
process  is  finite  and  terminates  with  a global  solution  of  the  approximate 
problem. 

3.1  Transshipment  on  Affine  Networks  - Economic  equilibria  on 
certain  affine,  multicommodity,  transshipment  networks  were  first  studied  by 
the  regional  economists  Taykayama  and  Judge  [26]  using  quadratic  programming. 
Recently,  Asmuth,  Eaves,  and  Peterson  [4]  have  constructed  a more  general 
approach  that  utilizes  the  economic  equilibrium  conditions  directly  without 
first  passing  to  a quadratic  programming  problem.  A brief  discussion  of 
their  model  and  the  solution  to  the  sample  problem  presented  in  their  paper 
follow. 

The  transshipment  problem  can  conveniently  be  represented  by  a directed 
graph  (N,L)  with  a finite  number  of  nodes  (members  of  N)  and  links  (members 
of  L)  on  which  a finite  number  of  commodities  can  be  transported.  Each  node  i 
in  N represents  the  set  of  producers  and/or  consumers  at  a specific  spatial 
location;  and  each  link  ij  in  L represents  a specific  facility  for  trans- 
porting commodities  from  some  node  i to  a different  node  j . (In  particular, 
we  assume  that  there  are  no  loops;  i.e.,  no  links  connecting  a given  node  to 
itself.)  Each  link  is  aligned  to  coincide  with  a direction  of  possible  trans- 
port; therefore,  at  least  two  links  must  connect  nodes  between  which  com- 
modities can  be  transported  in  either  direction. 
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The  nodes  are  enumerated  in  any  order,  consecutively  beginning  with 
one  - as  illustrated  by  the  graph  in  Figure  1. 


Figure  1.  Sample  Transshipment  Network. 


Let  p^  and  p^  be  n vectors  denoting  the  unit  demand  price  and 

the  unit  supply  price  of  n commodities  at  node  i , and  let  p1-*  e Rn 
denote  the  cost  of  transporting  each  of  these  n commodities  over  link  ij 
Affine  relations  are  assumed  between  prices  and  quantities;  i.e.. 


p1  = A^x*'  + a"*-  for  each  node  i £ N 


(5) 


where  p1  = (P^*Pg)  » A1  and 


are  given  constant  matrices  and  vectors 


(which  arose  in  [4]  from  inverting  the  difference  between  given  supply  and 
demand  quantities,  originally  expressed  as  affine  functions  of  price),  and 
x1  Is  an  n-dimensional  vector  representing  the  excess  quantity  of  each 
commodity  produced  by  node  i . It  is  also  assumed  that 
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p1^  = A1''x1'*  + a1'*  for  each  link  ij  c L 


where  A1"1  and  a1"1  are  given  constant  matrices  and  vectors  (which  arose 
from  describing  the  transport  prices  as  functions  of  transport  volumes),  and 

x1J  e Rn  denotes  the  quantity  of  n commodities  transported  over  link  ij  . 

The  quantities  are  constrained  by  the  nonnegativity  condition 

x1'1  > 0 for  each  link  ij  e L (7) 

and  the  commodity  conservation  condition 

x^  = T x*'*  - y xji  for  each  node  i £ N (8) 

j£N  j£N 

where  x1^  =0  if  ij  t L . Note  that  although  x1J  is  nonnegative  by 
virtue  of  the  choice  of  link  direction,  the  components  of  x*  might  be 
positive  or  negative  depending  on  whether  node  i is  a net  exporter  or 
net  importer  of  a particular  commodity. 


ship : 


The  price  stability  condition  on  p leads  to  the  following  relation- 


> p"l  for  each  link  ij  £ L 
*s  = d J 


To  see  this,  assume  that  (9)  is  violated  for  some  commodity  c . As  a 
consequence  some  economic  agent  would  find  it  profitable  to  purchase  as  much 
of  commodity  c as  possible  at  node  i and  transport  it  over  link  ij  for 
resale  at  node  j . This  would  clearly  be  an  economically  unstable  situation. 

The  final  relationship  needed  to  establish  equilibrium  is  the 
complementarity  condition 

< xij,  (p;j  - p^  - p1")  > = 0 for  each  link  ij  E L . (10) 

Cl  s 

This  condition  is  imposed  to  insure  that  no  positive  flow  x1^  will 
occur  on  a link  if  the  cost  p^  + of  a commodity  at  node  j exceeds  the 

price  p^  which  a consumer  is  willing  to  pay. 
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To  conform  with  the  notation  developed  in  Section  1.0,  a function 

g1-'  e Rn  equivalent  to  the  excess  demand  function  but  now  expressed  in 
terms  of  prices  rather  than  quantities  will  be  defined  by  the  following 
expression: 

g^  = p^  - p^  - p^  for  each  link  ij  e L 

Through  the  appropriate  substitutions  the  solution  (x,p)  to  the 
economic  equilibrium  conditions  (5)  - (10)  can  be  described  entirely  in 
terms  of  the  solution  (x,g)  to  the  linear  complementarity  conditions 


<x,g>  =0  g < o 

(li) 

0Q 

II 

1 

1 

< 

(12) 

where  x , g , and  v are  vectors  equal  in  size  to  the  number  of  links 
times  the  number  of  commodities,  and  M is  a square  matrix  of  comparable 
dimension  whose  components  are  given  in  Figure  2.  The  constant  v follows 
from  the  substitution  of  (5)  and  (6)  into  (10)  and  is  given  by 

v1J  = a"!  - alj  - a1  for  each  link  ij  e L 
d s 

For  purposes  of  illustration,  the  2-commodity,  5-link  network  shown 
in  Figure  1 has  been  considered  for  analysis.  The  reworked  data  for  this 
problem  is  displayed  in  Figure  3.  If  conditions  (11)  and  (12)  are  now  put 
into  the  format  of  (S)  we  get  a problem  of  the  form 
10 

min  £ (min(0,w.)  + x ) (A) 

x , w i=l  1 

w.  - M,x  + x.  =v. 
i i i l 

- M^x  < v^  i = 1,2 10 

x>0 

where  is  the  ith  row  of  M . 

The  algorithm  that  is  used  for  the  computations  does  not  solve  the 
original  problem  S,  but  constructs  an  approximate  problem  to  solve  by 
replacing  each  of  the  associated  functions  with  their  piecewise  linear 


* 
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convex  envelopes.  A related  problem  is  simultaneously  introduced  which  gives 
a sharper  underestimate  of  the  optimal  value  of  the  approximating  problem 
than  does  the  convex  envelope  problem.  It  is  this  related  problem  that  the 
branch  and  bound  procedure  solves  first  to  get  estimates  on  the  optimal  value 
of  the  approximating  problem,  and  to  set  up  new  problems  if  the  estimates  do 
not  yield  a global  solution. 

The  functions  defining  the  constraint  region  of  problem  A are  all 
linear,  and  hence  convex,  and  therefore,  will  not  be  replaced  in  the  approximate 
problem.  The  functions  associated  with  the  nonlinear  variables  w^  in  the 

objective  function  (i.e.,  min(0,w^),  i = l,...,n)  are  piecewise  linear,  but 
concave  and  will  be  replaced  in  the  approximate  problem  by  their  convex 
envelopes,  which  in  this  case  are  straight  lines.  This  is  illustrated  in 
Figure  4. 


Figure  4.  Convex  Envelope  of  min(0,w,.)  . 
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The  branch  and  bound  technique  proceeds  to  divide  the  domain  of  these 
functions  into  pieces  corresponding  to  their  linear  segments  and  separately 
solves  the  set  of  related  problems  in  which  the  nonlinear  variables  are  respec- 
tively limited.  When  every  function  is  piecewise  linear,  as  is  the  case  with 
(A),  we  get  an  exact  solution  to  the  original  problem. 

It  is  customary  with  branch  and  bound  methods  to  describe  the  algorithm 
in  terms  of  a branch  and  bound  tree.  The  nodes  of  the  tree  correspond  to  the 
related  linear  subproblems,  while  the  branches  of  the  tree  correspond  to  the 
set  on  which  the  branching  variables  are  defined.  A solution  is  obtained  when 
the  best  upper  bound  at  any  node  is  equal  to  zero.  That  is,  any  feasible 
point  yielding  a subproblem  value  of  zero  necessarily  satisfies  the  equilibrium 
conditions.  For  problem  A,  no  tree  developed  because  the  solution  was  obtained 
on  the  first  iteration  of  the  algorithm.  The  optimal  vector  x is  given  by: 

x1  = (0.2353,  0.7059)  x2  * (0.0,  2.2941)  x3  = (1.5294,  0.0) 
x4=  (1.0098,  0.2151)  x5  = (0.0,  0.0)  . 

The  formulation  of  (A)  required  the  addition  of  10  auxiliary  variables 
to  the  original  set  of  10  linear  variables.  The  former  were  each  divided  into 
two  intervals  for  the  purposes  of  branching.  This  division,  corresponding  to 
the  segments  of  the  piecewise  linear  functions  defined  for  these  variables, 
implied  that  any  branch  and  bound  tree  produced  by  the  algorithm  could  be  at 
most  10  branches  deep  and  that  no  variable  could  appear  more  than  once  along 
any  path.  In  theory,  it  might  have  been  necessary  to  solve  up  to  2^  - 1 
subproblems  before  reaching  a solution;  however,  the  fact  that  the  first 
subproblem  produced  an  equilibrium  point  underscores  the  computational  effi- 
ciencies that  result  from  having  available  at  the  outset  a means  of  independently 
checking  each  iteration  for  convergence. 

Each  subproblem  solved  by  MOGG  is  a linear  program.  When  the  excess 
demand  functions  are  affine,  the  solution  vector  necessarily  yields  a feasible 
point  to  the  original  problem.  The  upper  bound  associated  with  this  feasible 
point  will  always  be  greater  than  or  equal  to  zero,  but  generally  not  corres- 
pond to  an  equilibrium  solution.  Mangasarian  [18]  has  shown  that  the  linear 
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.. 


complementarity  problem  is  equivalent  to  a linear  program  whose  cost 
coefficients  are  dependent  upon  the  structure  of  M . The  similarity 
between  these  linear  programs  and  those  set  up  by  MOGG  when  M meets  cer- 
tain conditions  admits  the  possibility  that  MOGG  will  produce  an  equili- 
brium point  on  the  first  iteration.  Although  these  conditions  might  logi- 
cally arise  in  some  economic  contexts,  they  were  not  present  in  this  example 
and,  hence,  did  not  influence  the  rate  of  convergence. 

3.2  A Piecewise  Linear  Market  - This  example  [7]  provides  a simple 
explanation  of  how  a competitive  market  operates.  As  is  common  in  micro- 
economic  theory,  we  will  distinguish  among  individuals  according  to  the 
economic  functions  that  they  perform  or  on  the  basis  of  the  kinds  of  deci- 
sions they  make.  Thus,  a consumer  is  an  individual  (or  unit)  that  consumes 
commodities  and  supplies  inputs  to  production.  The  role  of  the  consumer  may 
be  defined  as  that  of  choosing  from  among  the  alternative  commodity  bundles 
available  to  him.  Similarly,  a producer  is  an  individual  (or  group)  that 
utilizes  inputs  to  produce  commodities.  The  role  of  the  producer  may  be 
characterized  as  that  of  choosing  from  among  the  alternative  input-output 
patterns  available  to  him.  The  same  individual  might  appear  in  the  economy 
both  as  a consumer  and  as  a producer.  Once  the  choices  are  made,  a state 
of  the  economy  is  defined. 

Under  certain  assumptions  (see.  Quirk  and  Saposnik  [22])  fora  market 
that  contains  n commodities,  m consumers,  and  £ producers,  the  aggregated 
(net)  amounts  of  commodities  demanded  and  supplied  for  any  vector  of  prices 
can  be  determined  by  a simple  summation  of  the  amounts  demanded  and  supplied 
by  individual  consumers  and  producers.  Thus,  given  the  price  vector  p , 
where  p = (P]_»P2»  • * * *Pn)  » we  can  wr*-te  x^  (p)  as  the  amount  of  the  ith 

commodity  consumed  (or  supplied  as  an  input  in  production)  by  the  jth 
individual  at  the  set  of  prices  given  by  p ; and  y^(P)  as  the  amount  of 

the  ith  commodity  produced  (or  used  up  as  an  input  in  production)  by  the  kth 
firm  at  the  set  of  prices  given  by  p . Then,  the  aggregate  (net)  consumption 
of  commodity  i by  consumers  is  given  by 
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xi(p) 


l 


j = l 


8iJ(p) 


1 )fi 


and  the  aggregate  (net)  production  of  commodity  i by  producers  is  given  by 

£ 

y, (p)  = l y,u(p)  i = 1,2, .. .n 
1 k=l 

We  tlien  define  x(p)  and  y(p)  as  point  to  point  mappings  from  Rn  into 
itself.  In  the  absence  of  any  initial  endowment  the  excess  demand  function 
can  be  written  as  g(p)  = x(p)  - y(p)  . 

The  sample  economy  under  consideration  contains  three  commodities, 
three  consumers,  and  three  producers.  The  associated  supply  and  demand 
functions  are  assumed  to  be  piecewise  linear,  and  are  given  in  graphic  form 
in  Figures  5 and  6.  To  conform  with  the  presentation  in  [7],  the  equilibrium 
quantity  rather  than  the  equilibrium  price  will  be  computed.  The  following 
notation  will  be  used: 

p = jth  producer's  supply  price  for  commodity  i 

Sij 

p^  = jth  consumer's  demand  price  for  commodity  i 

ij 

x^  = quantity  of  ith  commodity  consumed 


y^  = quantity  of  ith  commodity  produced 


where  i,j  = 1,2,3,  p = l Po 

1 J IJ 


x and  p , * } p , 

d . V d . . 

i j ij 


is  a function  of  the  consumption  variable 
is  a function  of  the  production  variable  y . 


An  equilibrium  point  will  exist  if  the  following  conditions  are 
satisfied : 


x > 0 , y > 0 

(13) 

x - y = 0 

(14) 

p , - p <0 
d rs  = 

(15) 

<x,  p , - p > = 0 
d s 

(16) 

15  - 


r 


T-386 


where  p^  and  p^  are  the  three-dimensional  market  supply  price  and  market 

demand  price  vectors.  The  first  condition  assures  feasibility:  the  second 
condition  clears  the  market;  the  third  condition  assures  price  consistency 
by  requiring  the  excess  demand  function  to  be  less  than  or  equal  to  zero; 
and  the  fourth  condition  is  Walras'  Law  and  reflects  the  following  circum- 
stances: if  , the  quantity  of  the  ith  commodity  being  purchased,  was 

positive  and  if  the  producers'  supply  price  p was  greater  than  the 

i 

consumers'  demand  price  p , then  the  producers  would  be  losing  money  and 

di 

begin  to  lower  y^  (=x^)  to  zero.  Such  a situation  would  be  economically 
unstable . 

Conditions  (13)  - (16)  can  now  be  written  as  a minimization  problem 
in  the  form  of  (P) . 

min{<-x  , (Pd_Ps)>^ 

p - p <0 

d s = 
x > 0 

In  order  to  recast  this  problem  in  the  form  of  (S) , the  consumption  and 
production  data  given  in  Figures  5 and  6 must  be  aggregated  over  their  respec- 
tive agents  to  obtain  the  market  demand  and  supply  curves  p,  and  p 

d s 

This  has  been  done  for  each  of  the  three  commodities. 


Commodity  1 


(-ZXj+18  ; xx  < l'l 

*1  (-0. 5x  +16. 5;  xL  > l| 


+ (6-3x2)  + (5-2X3) 


= Jllxl  * *1  2 

51  (7.5x3  - 3.5;  > if 


+ (2+x2)  + (2+2x3) 


18  - 
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Commodity  2 


p,  = (12-5x  ) + (14-2x  ) + (10-5x  ) 
d^  I Z J 


p = 4 + 

S2 


Commodity  3 


2 ’ 2 
<2+9;  x2 


p = 2 + (5-2x  ) + 
3 


'1  = 
;i  =• 


< n 

> + 3 

> ll 

-X-  + 8 ; 

x„  < l) 

3 

-2.5x2  +9.5; 

x > 17 

1) 

( 5x  ; x 

> + (3+x.)  + 

} 

l)  2 

(0.5x^  + 4.5;  x 

i( 


The  minimization  problem  in  its  separable  form  becomes 


min  } min  0,w.  + x.} 

. , i i 

x,w  1=1 

w.+p,  -p  +x.=0 
l *d.  *s.  l 

i l 

p - p <0 

d.  s.  = 

l l 


(B) 


i = 1,2,3 


x > 0 


where  p , and  p are  defined  above. 
di  Si 

Each  of  the  six  variables  in  (B)  is  nonlinear,  the  first  three 
(Wj,w2,w.j)  corresponding  to  the  auxiliary  variables  and  the  second  three 

(x^,x2,x,j)  to  the  original  problem  variables.  The  associated  functions 

are  piecewise  linear  and  contain  at  most  one  break  point.  This  means  that 
the  branch  and  bound  tree  can  be  at  most  six  nodes  deep  and  that  a maximum 

of  2^-1  subproblems  might  have  to  be  solved.  Once  again  though,  the 
algorithm  converged  on  the  first  iteration.  The  computed  best  upper  bound 
for  the  first  subproblem  is  zero  and  hence  the  solution. 
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If  the  algorithm  is  permitted  to  run  past  this  point  until  its  usual 
termination  conditions  are  met  a total  of  19  subproblems  * ill  be  solved. 
Figure  7 depicts  the  resultant  branch  and  bound  tree  which  serves  to  illus- 
trate both  the  advantage  of  knowing  the  optimal  value  at  the  outset  and  the 
amount  of  work  required  to  search  for  alternative  solutions.  The  known 
results  are  corroborated  at  node  9.2  where  roundoff  errors  have  produced  a 

-4  -5 

best  upper  bound  of  -0.2  x 10  and  a best  lower  bound  of  0.1  x 10  , 

a minor  contradiction.  The  two  numbers  adjacent  to  each  node  represent  the 
best  upper  and  lower  bounds  for  that  subproblem.  A bar  in  the  place  of  the 
best  upper  bound  indicates  that  no  corresponding  feasible  point  to  the 
approximate  problem  exists.  The  numbers  along  the  branches  refer  to  the 
branching  variables  associated  with  the  preceding  node,  and  the  + and  - 
signs  indicate  whether  the  particular  auxiliary  variable  was  permitted  to 
range  over  the  set  of  positive  real  numbers  or  negative  real  numbers,  respec- 
tively. The  bars  appearing  below  the  nodes  indicate  that  either  the  lower 
bounds  of  the  associated  subproblems  are  all  greater  than  the  current  best 
upper  bound  or  that  they  are  infeasible  and,  therefore,  cannot  contain  the 
solution. 

In  terms  of  the  actual  variables,  the  solution  vectors  are 
x*  = (2,1,1)  and  w*  = (-2, -1,-1)  . From  the  equality  constraints  in  (B), 
it  can  be  seen  that  x*  + w*  = 0 whenever  the  corresponding  excess  demand 
functions  are  binding.  By  tracing  the  convergent  path  backwards  from  node 
9.2  to  node  0.0  we  see  that  the  branches  that  fall  along  this  path  correspond 
to  the  nonpositive  orthant  of  w . The  subscripts  attached  to  the  branching 
variables  in  the  tree  denote  the  (closed)  intervals  over  which  the  original 
problem  variables  are  defined  for  all  subsequent  subproblems. 

3.3  A General  Market  - The  separable  programming  algorithm  works  by 
first  replacing  each  of  the  original  problem  functions  with  their  piecewise 
linear  convex  envelopes,  and  then  creating  a new  problem  to  solve  as  an 
ultimate  approximation.  From  this  approximate  problem  a series  of  convex 
subproblems  issue  that  are  set  up  and  solved  under  the  branch  and  bound 
philosophy.  If  the  original  functions  are  all  piecewise  linear  (but  not 
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Figure  7.  Branch  and  Bound  Tree  for  Piecewise  Linear 
Market  When  Forced  Past  Solution. 
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necessarily  convex)  then  solving  the  aggregate  of  subproblems  is  tantamount 
to  solving  the  original  problem  exactly.  Such  was  the  case  in  both  examples 
1 and  2.  In  this  example,  four  of  the  piecewise  linear  functions  in  Figures  5 
and  6 have  been  replaced  with  smooth  counterparts.  The  new  functions  were 
constructed  to  pass  through  the  points  (1,*)  and  (2,*),  and  are  given  by: 


16.516  e(-°*03175xl> 


"11 


33 


= -0.75  x^  - 0.25  x^  + 8 


= -3.75  x7  + 14.75  x, 
S11  1 1 


= 17.31234  log(x°'463+l) 
S22  2 


Substituting  these  functions  for  the  originals  in  problem  B leads  to  a new 
minimization  problem  that  can  be  written  as 


min  £ (min(0,w.)  + x.) 

w,x  i=l  1 


(C) 


subject  to 


, , „ 2 , ,,  , (-0.03175X, ) 

w^  + 3.75  x^  - 13.75x^  + 16.516  e 1 


4x2  - 4x3  = - 7 


w2  - 5x^  - 17.31234  log(x2‘4^3+l)  - x2 


5x3  = - 29 


; X1  = 1) 

j 

\ ~ 3x„  - < 

|0.75x2  + 4.25x3  - 8 ; x3  < l| 

i O M 

3;  Xj.  > 1) 

. / 

l0.75x3  - 0.25x3  - 3.5;  x3  > lj 

= - 4 


. .r  2 . , c1,  (-0.03175  x ) 

3.75  x^  - 14.75  x^  + 16.516  e 1 


4x2  - 4x3  < 


5x3  - 17.31234  log(x2'463+  1)  - 2x2  - 5x3  < - 29 


3x^-3;  x^  > 1 


- 3x2  - 


0.75x2  + 5.25x3  - 8 ; x < 1 

,2 


0.75x~  + 0.75x3  - 3.5;  x > 1 


< - 4 


x^  > 0 , x2  > 0 , x3  > 0 
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The  number  of  cuts  required  to  approximate  a piecewise  linear 
function  exactly  is  equal  to  the  number  of  segments  constituting  that  func- 
tion. When  the  function  is  smooth,  it  cannot  be  represented  exactly  by  a 
finite  number  of  linear  segments  but  can  be  approximated  with  arbitrary 
precision  by  increasing  the  number  of  cuts.  In  this  example,  six  cuts  were 
made  in  the  original  problem  variables  (x^.x^jX^)  over  the  closed  interval 

[0,3].  Owing  to  the  fact  that  the  cuts  were  evenly  spaced  every  half  integer, 
and  that  the  graphs  of  the  smooth  functions  pass  through  the  solution  points 
of  (B),  it  is  reasonable  to  expect  that  the  solution  to  (C)  would  coincide 
with  one  or  more  of  these  points.  This  indeed  was  the  case,  the  identical 
solution  x*  = (2,1,1)  , resulted  for  (C).  The  associated  branch  and  bound 
tree  is  shown  in  Figure  8.  The  algorithm  is  seen  to  have  converged  in  the 
tenth  stage  at  node  10.1  after  22  subproblems  had  been  solved.  This  con- 
trasts with  the  first  two  examples  where  the  solution  occurred  on  the  first 
iteration;  however,  in  each  of  these  three  problems,  the  first  feasible 
point  produced  by  the  algorithm  resulted  in  the  solution.  Finally,  we 
observe  from  Figure  8 that  at  the  tenth  stage,  the  best  upper  bound  and  the 
best  lower  bound  are  nominally  equal  implying  that  the  general  conditions 
for  optimality  have  been  satisfied,  so  the  algorithm  terminated.  If  an 
equilibrium  point  had  not  yet  been  found  at  this  stage,  it  would  have  been 
reasonable  to  conclude  that  none  existed  for  the  given  model.  The  other 
two  equilibrium  points  were  not  uncovered. 

3.4  Network  Traffic  Flow  - The  model  of  the  road  system  considered 
here  derives  from  the  notion  that  there  exists  a large  community  of  users, 
each  of  whom  takes  the  quickest  route  available,  given  the  actions  of  other 
users.  The  number  of  trips  taken  is  assumed  to  depend  on  the  time  required 
to  make  a trip,  while  the  travel  time  on  a particular  road  is  assumed  to 
depend  on  traffic  volume.  The  example  that  we  will  investigate  was  studied 
by  Asmuth  [3]  who  used  stationary  point  theory  in  conjunction  with  the  Eaves- 
Saigal  algorithm  [9]  to  obtain  a solution.  As  will  be  seen,  the  traffic 
flow  problem  closely  resembles  the  multicommodity  network  presented  in 
example  1. 
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To  formulate  the  model,  consider  a directed  network  (N,A)  with 
nodes  i in  N and  arcs  ij  in  A . For  each  arc  ij  , we  are  given  a 
delay  function  f which  expresses  travel  time  on  arc  ij  as  a function 

of  the  traffic  flows  on  the  arcs  of  the  network.  The  travel  time  along  arc 
ij  will  necessarily  depend  on  the  flow  on  that  arc,  but  may  as  well  depend 
on  flows  along  other  arcs.  For  example,  a two-way  street  could  be  modeled 
as  a pair  of  opposing  arcs  where  the  flow  on  one  of  the  arcs  might  affect 
the  travel  time  on  the  other. 

For  each  pair  of  distinct  nodes  i and  k we  are  also  given  a travel 

demand  function  g.  , which  expresses  demand  for  travel  from  i to  k as 
l , k. 

a function  of  travel  times  between  nodes  on  the  network.  Demand  for  travel 
from  i to  k will  depend  on  travel  time  from  i to  k as  well  as  on 
travel  times  between  other  pairs  of  nodes;  for  instance,  from  i to  some 
alternate  destination. 

Numerous  solution  procedures  have  been  proposed  for  computing  equilibrium 

traffic  flows  and  travel  times  on  the  network.  When  f is  integrable  and  g 

has  an  inverse  which  is  integrable  the  usual  approach  has  been  to  reformulate 

the  equilibrium  problem  as  a convex  programming  problem.  These  conditions 

will  be  met  if  each  f . depends  only  on  the  total  flow  along  arc  ij  , and 

il 

each  g.  , is  monotone  decreasing  and  depends  only  on  the  travel  time  from 
1 , k 

node  i to  node  k . Beckman  [6],  Florian  and  Nguyen  [14],  and  others  have 
addressed  tiiis  problem  under  comparable  conditions. 


In  practice,  the  demand  and  delay  functions  f and  g are  at  best 
empirical  fits  and  can  be  endowed  with  these  or  any  other  restrictions  which 
may  seem  useful.  Asmuth's  approach  does  not  depend  on  such  restrictions, 
but  only  requires  that  the  delay  functions  f^  be  positive  on  each  arc, 

that  the  travel  demand  functions  g.  , be  nonnegative  and  bounded  for  each 

i , k 

pair  of  nodes,  and  that  the  network  be  complete;  that  is,  a directed  path 
must  exist  from  every  node  to  every  other  node.  From  a strictly  analytic 
point  of  view,  we  will  require  only  that  it  be  possible  to  put  the  functions 
f and  g into  a separable  form.  However,  if  the  model  is  to  accurately 
reflect  the  properties  of  the  system  it  might  be  desirable  to  adopt  the  above 
restrictions. 
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The  mathematical  conditions  for  a user  equilibrium  are  presented 

below.  The  travel  time  from  node  i to  node  k will  be  written  as  t.  , 

i,k 

and  the  flow  on  arc  ij  with  destination  k will  be  written  as  y , . 

1 J » K- 


It  will  be  said  that  the  travel  time  vector  t 

equilibrium  if  the  following  conditions  hold: 

and 

flow  vector  y 

are  in 

gi,k(t)  = ? yij,k  " l yji,k 

i^k. 

i , kr.N 

(17) 

y > 0 

(18) 

Ci,k^  fij(y)  +tj,k 

Ck,k  = 0 

i^k, 

kf.N 

ijf-A, 

kcN 

(19) 

yij,k(fij(y)  + H'V*  = 0 

i^k. 

ijcA, 

krN 

(20) 

yij  £ yij,k 

ijr.A 

(21) 

Condition  (17)  is  the  conservation  of  flow  equation.  It  says  that  the  traffic 
leaving  node  i with  destination  k is  the  sum  of  the  traffic  arriving  at 
node  i with  destination  k and  the  traffic  originating  at  i with  des- 
tination k . Condition  (18)  says  that  traffic  flows  cannot  be  negative. 


Conditions  (19)  and  (20)  require  that  traffic  flow  by  the  fastest 
route  available.  In  (19)  we  require  that  t.  , not  exceed  the  minimum 

1 , K 

travel  time  from  i to  k given  the  flows  y on  the  network;  (20)  limits 

the  traffic  to  those  routes  which  achieve  this  minimum  travel  time.  Together 

(19)  and  (20)  imply  the  principle  of  minimum  travel  time.  This  says  that 

if  any  traffic  flows  from  i to  k ; that  is,  if  V y..  > 0 , then 

• 1 j * k 


t 


i,k 


min  (fij(y)  + tj>k) 


Equation  (21)  relates  the  basic  flows  to  the  total  arc  flows. 


It  may  be  useful  to  think  of  this  system  as  a multicommodity  network, 
where  all  of  the  traffic  destined  for  a particular  node  k is  a separate 
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commodity,  all  of  which  must  be  shipped  to  node  k via  the  network.  In 

this  way  g.  (t)  is  the  amount  of  commodity  k which  must  travel  from 
1 , K 

node  i to  node  k . This  trip  will  traverse  a path  of  arcs  from  i to  k 

Conditions  (17)  - (21)  can  be  put  in  the  form  of  problem  P as  follows 

min  I I < y.,  , . ('^(y)  " ^ k + ci  y)  > 
y>0  keN  ijeA  1J’k  1J  3,k  1>k 

t>0  i^k 

subject  to 


f.  - t.  . + t,  , < 0 

13  3>k  = 

i^k. 

ij  eA, 

i,k  " ^ yij ,k  + ^ yji,k  0 

i^k. 

i,keN 

ij  yij ,k  ^ 

ij  eA 

problem  in  the  form  of  (1.7) 

we  get 

min  l l {min(0,w  , ) + y ..  ,} 
y>0  k£N  ij eA  1J,k  1J,K 

t>0  itk 


subject  to 


w. . , - f . . - t,  . + t.  + y. . , = 0 i^k,  ijcA,  kcN 
ij,k  ij  j,k  i,k  13, k 

- f..  - t.  , + t.  . <0  i^k,  ij£A,  keN 

13  1 »k  i.k  = 

W,k-° 


yli  - £ yij,k  ■ 0 ij£A 

The  sample  problem  presented  below  is  from  [4]  and  is  based  on  the 
directed  network  shown  in  Figure  9. 
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Figure  9.  Sample  Traffic  Network 
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and  the  travel  demand  functions  are. 

80 

gl,2(t)  tL  2+l 


120 

gl,3(t)  t±  3+l 


g2,l^t^ 


C2,l+1 

100 

t2,l+^ 


if 

"2.1* 

C2,3 

if 

‘2.1  S 

t2,3 

I 80 


go  o (t ) 

^ , J 


h,3+1 


20 


C2,3+1 


g3,l(t)  ~ t 


60 


^3,2^^  ~ t 


3,1 

100 


+1 


3,2 


+1 


if 

‘2,1  = 

t2,3 

if 

C2 , 1 * 

t2,3 

When  more  than  one  function  value  is  given  at  a particular  point 
(e.g.,  ^2  1 = t2  3^  the  va^ue  8 is  the  convex  hull  of  the  two  values. 

In  this  case  some  of  the  travelers  from  node  2 will  go  to  either  1 or  3 
depending  on  which  is  closest.  If  the  travel  times  are  equal  then  those 
travelers  who  want  to  go  to  either  1 or  3 will  be  divided  between  the  two 
destinations. 


In  their  present  form,  the  demand  functions  g 


2,1 


and  g2  ^ exhibit 


an  implicit  dependency  on  the  travel  times  t2,l  ar)8  3 anc*  therefore, 

must  be  made  separable  before  the  equilibrium  problem  can  be  solved.  Although 
this  cannot  be  done  explicitly,  the  desired  result  can  be  achieved  by 
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considering  the  following  three  disjoint  partitions  of  t : 

C2,l  < t2,3  ’ t2,l  = C2,3  ’ t2,l  > t2,3  . 

The  mathematical  program  associated  with  each  of  these  partitions  comprises 
26  variables  and  27  constraints.  Of  the  26  variables,  12  are  of  the  type 
required  to  achieve  separability  of  the  functions  while  the  remainder  are 
defined  in  the  original  problem  statement. 

The  solution  was  uncovered  in  the  third  partition  at  the  8<4th  stage 
after  168  subproblems  had  been  solved,  and  once  again,  coincided  with  the 
first  feasible  point  found.  The  resultant  branch  and  bound  tree  is  not  dis- 

I 

played  because  of  its  extensive  length,  but  the  final  computations  are  high- 
lighted in  Table  1 along  with  the  results  obtained  by  the  Eaves-Saigal 
algorithm.  The  minor  discrepancies  observed  between  the  variable  and  func- 
tional  values  computed  by  either  method  can  be  attributed  to  the  grid  size 
superimposed  on  the  algorithm  and  are,  hence,  subject  to  control.  Finer 
resolution  is  strictly  a matter  of  increasing  the  number  of  grid  points  pre- 
scribed for  the  original  nonlinear  variables  and  solving  a proportionately 
larger  problem. 

4.0  Conclusions 

The  computation  of  equilibria  plays  a major  role  in  the  analysis  of 
economic  and  transport  systems.  Whenever  the  equilibrium  problem  can  be 
formulated  as  a set  of  complementarity  equations,  we  have  shown  for  those 
cases  where  the  original  functions  are  implicitly  separable,  that  nonconvex 
programming  can  be  used  to  obtain  a solution  to  either  problem.  A general 
algorithm  based  on  branch  and  bound  techniques  was  adapted  to  perform  the 
equilibrium  computation.  Unlike  the  usual  nonconvex  program  though,  where 
the  solution  is  recognized  only  when  equality  is  achieved  between  the  best 
upper  and  best  lower  bounds,  an  independent  check  can  be  made  for  the  solu- 
tion at  each  iteration  because  the  best  upper  bound  is  known  at  the  outset. 

As  our  computational  experience  demonstrates,  this  enhancement  markedly 
increases  the  efficiency  of  the  algorithm. 


\ 


However,  the  fact  that  a numerical  procedure  will  terminate  with 
the  correct  answer  in  a finite  number  of  iterations  is  no  guarantee  that 
it  will  be  of  any  practical  use.  The  combination  of  method  and  algorithm 
under  study  derives  its  tentative  usefulness  from  the  observation  that  for 
most  problems  investigated,  convergence  occurred  in  a far  smaller  number  of 
iterations  than  theoretically  possible.  The  results  have  been  especially 
encouraging  for  problems  of  larger  dimensions;  and  in  all  cases,  the  equilib- 
rium solution  coincided  with  the  first  feasible  point  found  by  HOGG. 

The  affine  equilibrium  problem  or  linear  complementarity  problem 
holds  a particular  interest  because  of  its  unique  structure  and  implicit 
relationship  to  an  equivalent  linear  program.  Because  of  the  similarity 
between  the  first  subproblem  set  up  by  MOGG  and  this  linear  program,  immediate 
solutions  are  often  obtainable  from  MOGG  at  little  extra  cost.  In  fact,  the 
additional  work  required  to  determine  the  equivalent  linear  program,  even 
for  relatively  small  problems,  is  often  more  burdensome  and  more  computa- 
tionally expensive  than  permitting  MOGG  to  run  beyond  its  first  subproblem 
to  a point  of  convergence.  A further  and  decided  advantage  of  MOGG  is  that 
it  will  solve  all  affine  equilibrium  problems  regardless  of  their  matrix 
structure.  By  contrast,  the  majority  of  alternative  procedures  available 
are  limited  in  their  application  to  a number  of  special  cases  which  do  not 
necessarily  arise  in  practice. 
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